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(N ■ Abstract. In this paper we report on a systematic search for a neutral boson in the mass range 

5^ \ between 5 and 15 MeV/c^ in the decay of highly excited nuclei. Its signature is found a deviation in 

■ the angular correlation of the e+e^ pairs from conventional internal pair conversion (IPC) resulting 

from of its two-body decay kinematics. With an e+e^ pair-spectrometer, a number of transitions 
. has been investigated in the a-nuclei ^Be, '^C and '^O, following Ught ion induced reactions at low 

bombarding energies, first at IKF in Frankfurt and during the last years at ATOMKI in Debrecen. 
\ Startlingly, in all isoscalar transitions excess e+e^ pairs are found at large angles with branching 

, ratios with respect to the total yield ranging from 10^^ to 10^^. If these deviations are all related 

to the two-body decay of an X-boson, this observation implies plural X-bosons. An analysis of all 
angular spectra with a boson search program, yields a pandemonium of more than ten candidate 
bosons. 



o 



o 



Introduction 



, Elusive neutral X-bosons can be emitted in a nuclear transition within the constraints 

^1 of spin-parity, isospin and energy-momentum conservation. Superallowed {L=0) emis- 

^ ■ sion is expected for a pseudoscalar particle in a 0^^0+ MO transition, a scalar 

particle 0+ in an 0+— s>0+ EO, an axial vector 1+ in an Ml and a vector I particle in an 
k>( ; El transition. Allowed {L=l) emission occurs for 0+ and 0^ particles in El and Ml re- 

■ spectively, and for 1+ and 1^ in El, Ml but also in E2 and M2 transitions. The signature 

^ of Z-boson emission is expected to be the two-body decay into e+e pairs superposed 

on conventional internal pair conversion (IPC) [1], but with a marked difference in the 
distribution of the e+e^ opening angle. 

Since some time a search is carried out for such a neutral particle in high-energy 
transitions in self-conjugate a-nuclei. Since 1990, experiments have been performed 
at the 'small' Van de Graaff accelerator of the Institut fiir Kemphysik in Frankfurt, 
Germany and since the year 2000 at the AVE cyclotron and the Van de Graaff accelerator 
at ATOMKI in Debrecen, Hungary. In all experiments the same compact spectrometer 
was used, consisting of eight E-AE detectors, allowing a good compromise between 
angular granularity, stopping power, range and pair- efficiency [2]. 

A signal which could originate from a boson with an invariant mass of 9 MeV/c^ has 
been observed in the Ml transitions depopulating the 17.64 MeV, 1+ level in ^Be to the 
ground- and first excited state. In Table 1 of refs. [6, 5, 7], the results are collected which 
are relevant to this anomaly for the a-nuclei ^He to ^°Ne, like the branching ratio Bx, the 
width Fx, the coupling strength ax and the nix for the different transitions. It appears 
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that isoscalar magnetic transitions and a possible 0^— >0+ transition exhibit an excess 
of e+e pairs at large opening angles, whereas no deviations occur in isovector El and 
Ml transitions. The deviations in isoscalar Ml transitions indicate an isoscalar character 
for such a boson with spin-parity 0~ or 1+, and coupling strength ax proportional to 
the isoscalar strength of the Ml transition. The sensitivity of these Z-boson searches — 
expressed in the ratio of Fx to F^varies between 10"^ and 10~^. 

In order to test this X-boson scenario, we started a dedicated search of e^e~ angu- 
lar correlations in isoscalar magnetic transitions with inelastic proton scattering experi- 
ments on ^^C and ^^O, the (/?, a) reaction on ^^F and the (c'He.pe^e^) reaction on ^^B 
and ^^N targets. The first aim was to check the signal observed [7] for the isoscalar Ml 
ground state transition from the T=0 level at 12.71 MeV in ^^C at a high level of signif- 
icance and the second aim to look for the 9 MeV/c^ boson in the 'forbidden' O^^O"*" 
(MO) transition at 10.96 MeV in ^^O. The mere observation of e+e^ decay to the ground 
state imposes a straightforward bound on an X -boson emission width. 



Experiments 

Experiments are carried out at ATOMKI, in Debrecen, Hungary at the cyclotron and at 
the Van de Graaff accelerator. The first four experiments (DEB1-DEB4) were dedicated 
to the ^'^C{p,p'e^e-y^C* and the the ^^0{p,p e^e-^^O* reaction at 16.5 MeV at the 
cyclotron using only the six 'small' detectors. For calibration the ^^F(p, ae+e^)^^0* 
reaction was used at 3.5 and 5.5 MeV. In the DEB5-run, the ''^¥\(^He,pe'^e Y^O* 
reaction was studied at 2.3 MeV together with the ^^F(/7, ae+e~)^^0* reaction at 3.5 
MeV using the complete eight detector set. The DEB7 experiment was performed in 
an 'asymmetric' arrangement with five telescopes (two large and three small ones) 
and an additional detector for proton detection. The ^^B(^//e,pe+e~)^^C* reaction 
was measured followed by the ^^'H(^He,pe^e^Y^O* reaction for comparison. These 
reactions were found to favor the population of unnatural parity states [8, 9]. In this 
experiment a thin proton detector with an 1 1 % efficiency was implemented downstream. 
In coincidence with the populating proton group, e^e~ pairs from the different levels 
could be distinguished. 

In comparison with the p-capture spectra obtained at IKF where the decay of only one 
level is selected by the choice of the bombarding energy, the (/?, a), i^He^p) and {p^p') 
reactions have several entrance channels in the final nuclei. The resulting e~^e~ sum- 
spectra are composite, but nevertheless different sum-peaks can be distinguished. We 
clearly observe deviations in the angular correlation spectra of transitions with energies 
below 9 MeV, the aimed boson mass, which energetically cannot decay by this particle. 
In the particle scenario with (mx^E-sin{(o/2)), these deviations have to be explained 
by the emission of different particles with lower mass. With the observation of several 
deviations, a number of new particles has to be postulated. This also opens the exciting 
possibility of multiple boson emission in one transition. 

Analysis of the presumably composite spectra with the number of bosons, their mass 
and strength as free parameter is not trivial and time-consuming. We have solved this 
problem in an indirect approach. We converted the angular spectrum in a 'quasi '-mass 
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spectrum by comparing the ratio of experimental data-points, divided by the simulated 
IPC distribution, with the ratio obtained by simulation of e.g. 10^ bosons divided by 
an equal number of IPC's. Normalization was obtained by requiring an equal content of 
simulated to experimental X/IPC events by multiplying the theoretical ratio with a factor 
u, the transmission or inverse efficiency for boson emission. Boson masses were fitted 
from threshold to the transition energy in e.g 64 channel bins using a procedure. The 
inverse of the was taken to be the transition probability: P = {x^ / (iV — 1)) \ where 
N is the number of detector combinations. The value P can then be plotted as a function 
of assumed m,-, testing the likelyhood of a particle with the given mass being present in 
the spectrum. An upper limit for the branching ratio is obtained by multiplying P with 
the 'transmission' u. For three detector geometries 10^ or 10^ events have been generated 
for some 40 to 70 boson masses, with and without energy disparity (j) cuts, for transition 
energies of 6.05, 8.87, 10.96, 12.71 14.64, 17.23 and 17.64 MeV. Also for the IPC 
(EO, El, E2, Ml and M2) distributions, generally up to 10^ events were simulated with 
and without y-cuts. It was found that a proportional scaling between boson masses and 
transition energies could be applied, except for the low masses, where the transition from 
backward decay to forward two-body decay at /3cm = ^lab is not scaled correspondingly. 
A detailed description of the program will be published elsewhere [13]. 



The IKF experiments revisited 

To test the validity of the program we have first applied it to the earlier results 
at IKF. Spectra from the '^Li(/7,e+e~)^Be* reaction at 0.44 and 1.03 MeV and the 
^^B(p,e"'"e~)^^C* reaction at 1.6 MeV proton energy were considered. All 28 correlation 
angles were used in the fits. The mass distributions perfectly confirm the 9 MeV/c^ 
boson hypothesis with the peak at 9.23 MeV in Fig. Ic. The second peak at 10.95 MeV 
corresponds with the emission of the 9 MeV/c boson in the 14.64 MeV transition to 
the first excited state, being at 9.09 MeV/c^ in the 14.64 MeV laboratory system. No 
significant peak is seen at 9 MeV/c^ in Fig. la. The P- value for ^Be at 9 MeV/c^ is 
indeed a factor five larger than for ^^C as observed in the angular distributions. However, 
we also see peaks at the higher masses of 12.50 and 14.23 MeV/c-^ masses. 

A different picture is found for the isoscalar transition at 18.15 MeV. Here, the 
dominant peaks are at 12.42 and 14.55 MeV/c^. They are a factor 20 stronger than 
in the iso vector Ml, whereas the strength of the 9.23 MeV/c^ peak is about the same. 
The masses in Fig. 1 are indicated using the ground state transitions as two-body decay 
frame. Applying the transformation to 14.64, the masses become 2.99, 7.55, 9.26, 10.52 
and 11.86 respectively and in Fig. If (with 15.15) of 3.21, 4.78, 8.22, 10.37 and 12.15 
MeV. The emission of the same 9 MeV/c^ boson in the 17.64 (9.10) and in the 14.64 
MeV (9.23) was already suggested in [6] where approximately the same branching was 
deduced for this boson. Figs, lb. Id and If show the branching ratios Bx for boson to 7- 
ray emission. Integration over 9 MeV/c^ yields respective values of 1.9-10^^, 1.2-10^'^ 
and 4.0-10"4 within la in accord with the values of <2.3-10^^ (1.1±0.3)-10 and 
(5.8±2.2)-10~^ given in [3, 6, 5, 7]. The mass distributions for the transitions after 
applying symmetry cuts to the spectra, show dramatically larger variations than in Fig. 1 
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e^e" decay of the El and Ml transitions 
from the 17.23 1' in ^^C, 17.64 1*, T=l and 18.15 1*, T=0 in *Be 
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FIGURE 1. Figs, la, Ic and le show the probability distributions of bosons in the decay of the 
17,23 1", the 17.64 1+ and the the 18.15 MeV 1+ levels following the ^^B(p,e+e^y^C* reaction 
at 1.6 MeV and the ^Li(p,e^e^)^^0* reactions at 0.44 respectively 1.03 MeV as measured 
at IKF in 1992 and 1998 [3, 7]. Figs, lb. Id and le show the corresponding Bx-values i.e. 
branching ratios for boson to y-ray emission. The dashed lines in Figs. Ic and Id correspond 
with the El distributions from Figs, la and lb. 



at 9.0, 12.5 and 14.5 MeV/c^, in agreement with the expectations for two-body decay 
where the symmetrical decay is enhanced with respect to the primarily asymmetrical 
decay of Ml-IPC at large angles. Integration at 9 MeV/c^ yields Bx values of about half 
of those obtained from Fig. 1, in agreement with the expectations for two-body decay. 

It appears that the boson search program adequately describes the IKF-results by 
reproducing the previously found signal at 9 MeV/c^. However, surprisingly, it also finds 
signals for particles at the masses of 12.5 MeV/c^ and 14.5 MeV/c^. In the calibration 
runs: ^^F(;?, ae+e^)^^0* at 0.72 and 0.84 MeV just below and above threshold for the 
8.87 MeV 2^ level, two boson candidates are observed at 5.0 and 5.38 MeV/c^ with 
branching ratios as large as several percent. 
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Back to ATOMKI and further 



Because of the apparent feasibility of transformation of the angular spectra to mass- 
spectra and the unbiased character of the program we have systematically converted 
all available angular spectra from IKF and from ATOMKI and other sources into mass 
spectra. The oldest investigated spectrum was the celebrated angular correlation mea- 
surement of the EO transition at 6.05 MeV in ^^O dated back to 1949 [10] and the most 
recent one: the study of the e+e^ decay of the giant resonance decay at 22.23 MeV and 
the 15. 1 1 MeV isovector Ml decay in at Stony Brook [1 1] in 2003. 



The 4.44, 6.05, 7.66, 8.87, 10.96 and 12.71 MeV transitions 

In this paragraph, we discuss the mass-distributions of six different transitions in ^^C 
and 1^0, selected from results of the ^^C{p,p'e+e-y^C* and the ^^0{p,p'e+e-y^O* 
reactions at 16.5 MeV and from the ^^F(p, ae'^e~)^^0* reaction used for calibration. 

The mass-distributions for symmetrical (S) pairs are shown in Fig. 2, with the P- 
distribution in the left column and the branching ratio in the right. In all shown data, 
disparity cuts of j<0.5 have been applied. For every transition more than one spectrum 
is available showing in essence the same features as Fig. 2. In particular, for the 6.05 
MeV EO transition, data from twelve experiments have been analyzed. Eight of them 
show strong peaks at 5.0 and 5.3 MeV/c^. The four spectra where they are absent, stem 
from experiments where the maximum correlation angle is smaller than 1 10*^. 

The data are discussed from bottom to top. right to left. In the 12.71 MeV transition, 
we clearly recognize the boson at 9 MeV/c^, the original goal for the experiments in 
Debrecen. The branching ratio is about 8-10^"^ in agreement with the results of [6, 5, 7]. 
However, there are two new candidates at 10.39 and 11.19 MeV/c^, which are stronger 
than the one at 9 MeV/c^. The Ml data are compared (bottom left) with a measurement 
of the angular distribution of the forbidden 0+-^0+ MO transition at 10.96 MeV in ^^O, 
where a strong signal is seen at 7.85 MeV/c^. The X/IPC value is as large as 0.8 and 
for the Cy<l) cut (not shown here) it is even 2, suggesting that the decay via the 7.85 
MeV/c^ boson is the dominant decay mode in this MO transition. The two peaks at 8.99 
and 9.84 MeV can be ascribed to reflection from bosons of 7.35 and 7.85 MeV/c^ both 
occurring in the decay of the neighboring 8.87 MeV/c^ which could experimentally not 
be completely resolved from the 10.96 MeV sum-peak. It seems appropriate to assign, 
from the superallowed decay, the 7.85 MeV/c^ as a pseudoscalar particle and for the 
same reason the signals at 9.23 and 10.39 MeV/c^ as axial vector particles. A vector 
character would be preferred for the 11.19 MeV/c^ signal. 

In the 8.87 MeV M2 transition we see a boson candidate at 7.35 with the spin 1 and 
possibly the 7.85 0~. The 6.32 MeV/c^ peak could be due to reflection of the 7.85 in the 
10.96 transition. However, the signals at 2.25. 3.38 and 5.04 MeV/c^ in this transition 
are all viable boson candidates. Also in the mass-distributions of the two EO transitions 
at 6.05 and at 7.65 MeV, multiple resonances are displayed. A candidate boson at 3.36 
MeV/c^ is seen in all spectra. The occurrence in the EO, E2 and M2 transitions (allowed 
decay) could only be compatible with the spin-parity of 1~ for this particle. With the 
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FIGURE 2. From top to bottom (left) are shown the probability distributions of symmetrical (S) 
pairs in the 4.44 (E2), 6.05 (EO), 7.65 (EO) and 8.87 (M2), 10.96 (MO) and 12.71 (Ml ) transitions 
in ^^C and ^^O. Right are shown the corresponding branching ratios for boson emission with 
respect to EO-IPC (X/IPC) for the 6.05 MeV EO and the 10.96 MeV (MO) transitions and to 
y-raysfor the E2 at 4.44 MeV, the M2 at 8.87 and the Ml at 12.71 MeV. The Bx-value for the 
7.65 MeV EO is shown with respect to the parallel 0+— >2+ 3.20 MeV E2-y-ray. 



simultaneous presence of the candidates below 4 MeV in the EO and El transitions the 
spin is restricted to 1^. Likewise, the resonance at 5 MeV/c^ must, from its presence 
in the EO's and the M2 transition, have spin-parity 1 . With its appearance in both EO 
transitions and its absence in the M2, the 4.48 MeV/c^ candidate must have a scalar 
0+ character since a 1^ particle would also have shown up in the M2 transition. With 
some optimism, this figure indicates the existence of two scalar X-bosons at 4.5 and 5.3 
MeV/c^ and six vector bosons at 2.5, 3.3, 4.0, 5.0, 6.3 and 7.3 MeV/c^. 



Anomalous Internal Pair Conversion 



Signaling Elusive Light Neutral Particles 



6 



Conclusion and the future 



It appears that in every isoscalar transition considered with energies from 4.44 to 
12.71 MeV deviations from IPC are observed. This anomaly is explained in an isoscalar 
particle scenario which would require a multitude of bosons. The finding could have the 
bizarre consequence of a new type of light neutral matter instead of or/and the discovery 
of a new gauge boson. 

A second generation of dedicated experiments is required to test this exciting pos- 
sibility. The properties of every boson candidate have to be checked in great detail. A 
first experiment in that direction [13] has been undertaken in the DEB7 run, where the 
two-body decay properties of the bosons were verified, by mounting the central detector 
combination at an asymmetrical position. An alternative approach is being performed 
[14] at ATOMKI to study the angular correlations with an order of magnitude higher 
angularity, using wire chambers in front of the dE / dx-detectors. 



This paper has been presented at the International Symposium on Exotic Nuclear 
Systems, ENS'05, at ATOMKI, Debrecen, Hungary, 20-25 June 2005. 

I am indebted to my colleagues: J. van Klinken (KVI, Groningen), H. Bokemeyer 
(GSI, Darmstadt), K. Bethge, O. FrohUch, K.A. MuUer, K.E. Stiebing (IKF, Frankfurt), 
R. van Dantzig, T.J.Ketel (NIKHEF, Amsterdam), K. Griffioen (William and Mary, 
Williamsburg), A. van der Schaaf (Universitat Ziirich, Ziirich), M. Csatlos, L. Csige, 
Z. Gaczi, J. Gulyas, M. Hunyadi, A. Vitez, A. Krasznahorkay Jr. and A. Krasznahorkay 
(ATOMKI, Debrecen). 
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